Mapping local Hamiltonians of fermions to local Hamiltonians of spins 
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We show how to map local fermionic problems onto local spin problems on a lattice in any 
dimension. The main idea is to introduce auxiliary degrees of freedom, represented by Majorana 
fermions, which allow us to extend the Jordan-Wigner transformation to dimensions higher than 
one. We also discuss the implications of our results in the numerical investigation of fermionic 
systems. 

Many of the interesting strong correlation quantum phenomena that appear in Nature can be accounted for in 
terms of discretized models, i.e., in lattices. The corresponding effective theories often deal either with fermions, 
which can hop from one lattice site to another one as well as interact among them, or spins, which sit at different 
lattice sites. For local Hamiltonians in one spatial dimension, the Jordan-Wigner transformation (JWT) jj] allows 
us to map a local theory of fermions onto a local theory of spins, and thus those systems exhibit related physical 
phenomena. In higher dimensions, however, this transformation gives rise to non-local interactions between spins, 
and thus fermions and spins are expected to exhibit very different behavior, at least when they are described in terms 
of local Hamiltonians. 

In this paper we will show how one can always exactly map a local theory of fermions to the low energy sector of a 
local theory of spins in a lattice of arbitrary dimension. To achieve this, we introduce extra degrees of freedom in the 
form of Majorana fermions, which interact locally with the original ones, and which allow us to perform a JWT in 
higher dimensions than one. Similar ideas were previously used by Bravyi and Kitaev in the context of simulating 
fermionic quantum-computational circuits with qubits; there they showed how local fermionic gates can be simulated 
with local gates on qubits if the quantum states are rescricted to live in a specific sector of the Hilbert space that is 
determined by nonlocal gauge conditions. Also, R. Ball 3] has shown how these gauge conditions can be transformed 
into local terms of the Hamiltonian, except at the boundary where nonlocal terms remain. We present a mapping 
which has the advantage that the low energy sector is completely determined by local terms in the spin Hamiltonian. 
The spin Hamiltonians obtained are exactly of the type studied by Wen et al. where the opposite question was 
studied, i.e. how fermionic behavior could originate out of the low energy sector of a bosonic Hamiltonian. 

The question analyzed in this paper is clearly of fundamental interest. However, it may also have important 
implications in the numerical investigations of fermionic lattice problems, e.g. the Fermi-Hubbard model in two 
dimensions. In fact, our original motivation was to develop a method to represent fermionic problems in terms of 
local spin problems, so that we could use the algorithms we have recently developed to deal with those systems 
in dimensions larger than one [f|. In that case, a direct numerical simulation in terms of the fermionic degrees of 
freedom was not a good choice, since in Fock space one always obtains non-local terms (which are directly related to 
the Jordan-Wigner transformation), and ground states of generic systems with long-range interactions do not have 
to obey e.g. the area law 1 and other characteristic properties that allow for efficient simulations: ground states of 
local Hamiltonians are very special in the sense that they are completely determined by their local reduced density 
operators, and hence a variational method will work fine if it allows to approximate local properties well 0. This 
is one of the reasons why it is important that the low energy sector of the spin Hamiltonian obtained is determined 
by local terms, and does not contain nonlocal gauge conditions. Apart from that, our results may also prove useful 
for other standard simulation methods, such as quantum Monte Carlo, since there one also effectively simulates spin 
Hamiltonians instead of fermionic ones directly. 

From the mathematical point of view, the natural formalism to describe a state of fermions is second quantization. 
Within this language, the most general state of N fermions can be parameterized as 
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where a\ are creation operators obeying the usual anticommutation relations, and where |f2) denotes the vacuum. 
Note that, from the moment we consider second quantization, we have to enumerate the fermions and remember this 
ordering. After this choice, the associated Hilbert space (or Fock space) in which we can describe quantum states 
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becomes a tensor product of spin 1/2 's or qubits by the identification 

Ni,«2, = (4) (4) ■••( a ]v) 

So it is indeed simple to represent a collection of N fermions in a Hilbert space of N spins || . When trying to simulate 
fermionic systems, one has to work with the vector Ci lf ...j, N which is defined on this Hilbert space of spins. 

The obvious question is now of course how fermionic Hamiltonians and creation/annihilation operators look like in 
this spin representation. Let us e.g. consider the representation of the annihilation operator a k : 
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This can readily be derived by using the anticommutation relations. The effective Hamiltonian described in spin- 
space can now easily be calculated. The transformation we did is exactly equivalent to the JWT which maps com- 
muting spin operators to anticommuting fermionic creation/annihilation operators. From this point of view, the 
Jordan- Wigner transformation is not just a very useful trick but a natural consequence of the formalism of second 
quantization. The Hamiltonian as defined on the Hilbert space of spins is thus obtained by carrying out a JWT on 
the fermionic Hamiltonian. 

Let us illustrate this transformation for some specific example. Due to supcrsclcction rules, physical observables 
always contain an even number of creation/annihilation operators, and thus in the following we will consider fermionic 
observables of this sort. Here we recall how typical operators transform (we neglected constant terms): 

a{a k -» a% 
a\a k+n + c.c. -» of Sfe,n0f +n + a y k S kin u v k+n 

where S k _ n — ®i= k +i a 'i ^ s ^ ne string operator, which is the most distinctive feature of this transformation and arises 
when considering hopping terms between distant fermions. 

The central question of this paper is now the following: let us assume we have a system of fermions on some lattice 
with a specific geometry; is it always possible to map a local Hamiltonian of fermions to a local Hamiltonian of spins 
in this way? The answer to this question could clearly depend on the geometry and the type of interactions. Let us 
first consider the 1-dimensional case; one immediately sees that one can always choose the ordering of the fermions 
such that locality is preserved by the JWT. The situation seems to be very different in 2 dimensions. Let us, e.g., 
consider a translational invariant model with nearest neighbor interactions and hopping terms on a square lattice. 
Due to the string of a z operators, there exists no ordering of the fermions such that all local hopping terms in both 
horizontal and vertical direction map to local operators in the spin representation. The same problem arises in all 
geometries different from a 1-dimensional chain. In the following we will show how locality can nevertheless always 
be achieved by introducing extra degrees of freedom in the form of extra fermions. 

For simplicity, let us consider a 2-dimensional square lattice of spinless fermions with N columns and open boundary 
conditions. We enumerate the fermions as in Fig. 1 and consider interaction terms of the form 

a\a 3 (1) 
afaioja.,- (2) 

where either i = j, either i and j label two geometrically neighboring fermions. The JWT will map all operators of 
the form (|SJ) to local spin operators, but the hopping terms of the form Q will only result in local spin terms when 
j = i ± 1 (i.e. all horizontal hopping terms and halve of the vertical ones at the edges). So there is clearly only a 
problem with the vertical hopping terms. 

The idea is now to introduce an extra fermion bi for each mode fij, and impose that these extra fermions are in the 
ground state |x) of the Hamiltonian 

n avx = J2 Pa = E ~& + ~ 
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Figure 1: Enumerating fermions on a square 2-D lattice. 



(a) (b) 
Figure 2: Two possible graphs specifying different Hamiltonians TCaux for the auxiliary fermions. 



where the set {i,j} includes only pairs that correspond to the directed edges in Fig. 2(a) from vertex i to j 
(l4|. It can readily be verified that all the commute, P?. = 1, and that the gap of the Hamiltonian TL aU x is 2. Its 
unique ground state has the property that 



V(i, j) S {i,j}, 



Pn\x) = \x)- 



Actually, we can define the operators Cj = (6j + b\) and di — —i(bi — SJ) as two Majorana fermions making up the 
usual fermions Q- 

Next we consider the system consisting of the original fermions hi and of the auxiliary system together; geometrically, 
we assume that each site is occupied by two fermions, the original hi and the new one bi>. We leave the Hamiltonian 
on the auxiliary modes Tiaux unchanged, but alter the vertical hopping terms of the physical Hamiltonian in the 
following way: 



al&j -> a\ajP vr = 



Clearly, all these vertical hopping terms commute with H. aU x, and hence the Pij are constants of motion (i.e. +1), and 
as such the ground state of the full system will be the tensor product of the original ground state and |x) 0. But the 
advantage of the new formulation is immediately clear when the fermions are reordered in the way 1, 1', 2, 2©.. when 
doing second quantization. As a result, the hopping terms in the vertical direction will now correspond to quartic 
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interaction terms which will remain local by doing a JW transformation: the string operator of k' will cancel the one 
of k. Note that the configuration of the auxiliary fermions shown in Fig. 2(a) is only one possible choice and many 
other loop structures would lead to similar results (one possibility would e.g. be to have only loops with 2 vertical 
edges such as in Fig 2(b)). 

But the problem is not quite solved yet, because the JW-transformation will turn TL aux into a nonlocal Hamiltonian. 
Note however that the JW does not map the vertical hopping at the alternating edges of the lattice hopping into 
nonlocal ones. Consider e.g. the pair (N, N + 1) and associated Pn,n+i at such an edge (here N is the number of 
columns; N = 4 in Fig. 1); let us now make the following substitution in Ti. aux : 

Pn-i,n+2 — ► Pn-i,n+2Pn,n+i 

Pn~2,N+3 — ¥ Pn-2,N+3Pn-1,N+2 
P\,2N — * Pl,2NP2,2N-l 

and similar for all other rows. Due to the fact that all P commute, this substitution does not affect the ground state 
and the gap of Ji a ux- The advantage of using this new Hamiltonian is however obvious: the JW-transformation will 
map it to a local Hamiltonian with local plaquette terms: the string of cr z -operators generated by cv will be cancelled 
by the ones of c/j+jy. So we indeed managed to transform a local Hamiltonian of fermions into a local Hamiltonian 
of spins! One could argue that this indicates that fermionic systems potentially represent the low-energy sector of 
bosonic systems (see e.g. Q). 

Let us illustrate the kind of terms that will appear in the spin Hamiltonian for the case of the square N x N lattice 
of spinless fermions. Without loss of generality, we assume N to be even. We arrange the auxiliary fermions as in 
Fig. 2(a), and imagine that we have two superimposed layers of N x N qubits, the first layer corresponding to the 
physical qubits and the second one to the auxilliary ones. We label the Pauli spin operators acting on the first layer 
as Xk i, Yk.i, Zk,i and on the second layer X k ,i, Y k j, Z k ,i with 1 < k < N labelling the column and I the row. Let us 
systematically show how all terms in a fermionic will look like in the spin picture. 

First of all we observe that the following terms transform trivially: 

("l' ,i'®k> ,l'J —> (Zkd - 1) {Z k >,v - 1) 

Horizontal hopping terms transform as 

+ ^k+i,J® fc > ! — * (Xk,lXk+i,i + YkjYk+i.i) Zk,i 

Vertical hopping terms become 

a\ tl a kt i + -i +a\ tl+1 a k ,i -> (X k jX k j +1 + Y k jY k j + i) (^X kt iY k j +1 ^ (k = odd) 

-» (X Kl X kA+1 + Y kA Y Kl+1 ) (Y ktl X kt i +1 ) (k = even) 



depending on whether k is even or odd. 

The gauge conditions can be imposed by including the following 6-body local terms in the spin Hamiltonian: 



^ (Zk+i,iZk,i+i) \ Yk,lYk+\,lYk,i+iYk+\,i+x 

1— odd, k— odd 

^ {Zk+l,lZk,l+l) (Xk,lXk+l,lXkd+lXk+\d+l 



1— odd, k— even 



^ {Zk,iZk+\,i+x) ( Yk,iYk+\jYkj+iYk+\,i+x 

1— even, k— odd 

^ (Z k ,lZk+lJ+l) (Xk,lXk+l,lXk,l+lXk+l,l+lj 



1— even, k— even 
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It can easily be checked that these terms commute and also commute with the hopping terms. Finally, boundary 
conditions are required to fix the gauge uniquely (note that an ordering of the fermions and the corresponding JWT 
was crucial to make these terms local): 



k=odd 
k=odd 

E ^N,k^N,k+lXN,k+l 
=odd 

l^l,fc+l 



The sum of all those terms (of course with the appropriate prefactors of the original model) is the local spin 
Hamiltonian whose low energy sector is equivalent to the original fermionic model. Let's try to get some insight into 
all these terms. The hopping terms give rise to the well known XX + YY interactions, but these are mediated by a 
Z term in the horizontal case and a XY term in the vertical case; it is intriguing how the noncommuting nature of 
these extra terms should be related to the topological features of fermionic systems. Let us now have a closer look at 
the gauge conditions. Let us first assume that all spins in the first layer are completely polarized in the Z-direction; 
this would correspond to a completely trivial Hamiltonian for the original fermions. The qubits in the second layer 
would then interact via plaquette interactions of the form XXXX or YYYY, as all ZZ expectation values are equal 
to 1 . It can readily be checked that the effective Hamiltonian acting on the second layer is then exactly equal to the 
celebrated toric code Hamiltonian of Kitaev [i| if one rotates the lattice over 7r/4: the topological features of the toric 
code states seem to be exactly the ones needed to mediate the desired mapping from local interactions of fermions 
to spins. It would be very interesting to investigate this connection between the topological features of toric code 
states and emergent fermions further, and to arrive at a more intuitive explanation for the proposed mathematical 
construction. 

This construction can readily be generalized to fermions with spin and any regular lattice geometry; in e.g. a 3-D 
cubic lattice, one would need two extra fermions per site (i.e. 4 Majorana fermions). The Hamiltonian characterizing 
the auxiliary fermions can in this 3-D case still be made local by multiplying four commuting Cjdj-terms with each 
other. A triangular lattice can e.g. be obtained by embedding it into a square lattice; in this particular case some 
fermions would not be coupled to any other ones, and extra Majorana fermions guiding extra horizontal hopping terms 
have to be introduced. One can also include long-range hopping terms by multiplying them with multiple factors Py. 
A 2-D square lattice of fermions with an additional spin 1/2 degree of freedom could either be treated by considering 
two layers of the spinless system, or by working with spin 3/2 objects instead of qubits. 

An interesting problem arises when the same analysis is carried out on systems with periodic boundary conditions. 
Let us e.g. consider the case of a 1-D ring; here we can add two Majorana fermions to each original fermion, and use 
the link between the first unpaired Majorana fermion and the last one as a tool to make the original interaction term 
between the first and last fermion quartic. However, the JW-transformation will make the interaction term between 
these two Majorana fermions nonlocal, and there does not seem to be a way around it: the topology of the ring enforces 
us to have one nonlocal term in the final spin Hamiltonian jlq . Similarly, in the case of a cylinder one nonlocal term 
will appear, and in the case of a torus two nonlocal terms. In other words, we can speculate that the ground state of 
the corresponding spin system has the property that it cannot be the unique ground state of a local spin Hamiltonian: 
the associated local Hamiltonian (i.e. without nonlocal term) will possibly have a two-fold degenerate ground state, 
and a nonlocal operator which is topologically nontrivial will be needed to specify the sector. This phenomenon has 
indeed been identified in the case of the toric code states |S( of Kitaev, and is related to topological quantum order 

Instead of regular lattices, one can also ask the question whether a similar mapping is possible in the case of a star- 
like geometry which e.g. occurs in multichannel Kondo-related impurity problems [Illll2j . Let us therefore consider 
an impurity coupled to two different kind of fermions each having an extra spin degree of freedom. Enumerating the 
different fermions as in Fig. 3 (here fermions nr. 6 and 7 correspond to the impurity), one finds that again a long range 
string operator will appear between the second type of fermions and the impurity (i.e. between {6, 17} and {7, 18}). 
As in the previous case, we introduce two pairs of Majorana fermions and multiply the problematic hopping terms 
with the appropriate Majorana operators. This will indeed make the hopping terms local after the JWT. The problem 
in this case however is that after the JW there is no way to impose the identities ic^dk = +1 for both pairs in a local 
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Figure 3: Enumerating fermions in the case of a 2-channel Kondo-like impurity problem; 1-5 represent the spin up fermions of 
the first channel, 6-7 the impurity, 8-12 the spin down fermions of the first channel, and 13-17 and 18-22 respectively the spin 
up and down fermions of the second channel. 



way. But here one can easily see that there is no need to impose that constraint: the effective Hamiltonian after the 
JW is block diagonal, one sector corresponding to ickdk = +1 and the other one to ickdk = — 1; the second block is 
equivalent to the first one under local unitary operations. Indeed, the local transformation a, — > exp(i7ra|ai)ai on all 
fermions of one type will leave the interaction terms between them invariant, and will only lead to a negative hopping 
term when coupling them to the impurity. So the ground states in the two sectors are related by local unitaries: 
no gauge fixing interaction term needs to be imposed, and when doing simulations on the spin Hamiltonian one can 
always check the sector one ends up in. So it is again enough to diagonalize a local spin Hamiltonian. 

In conclusion, we have shown how it is possible to map a local problem of fermions onto a local problem of spins 
by using auxiliary Majorana fermions. This leads to a 2-dimensional version of the Jordan Wigner transformation. 
We believe that, apart from its fundamental interest, this map may become very useful in the context of numerical 
description of fermionic systems, since in that case local spin Hamiltonians seem to be better approximated. An open 
question is whether one can use this mapping to solve certain interesting spin Hamiltonians exactly by transforming 
it back into free fermions, as it is done in 1-dimensional investigations and in the recent work of Kitaev |l3| . 
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